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ABSTRACT 


We  extend  the  1ST  for  the  Benjamin-Ono  (BO)  equation,  given  in  Ref.[l],  in  two 
important  ways.  First,  we  restrict  the  1ST  to  purely  real  potentials,  in  which  case, 
the  scattering  data  and  the  inverse  scattering  equations  simplify.  Second,  we  also 
extend  the  analysis  of  the  asymptotics  of  the  Jost  functions  and  the  scattering  data 
to  include  the  nongeneric  classes  of  potentials,  which  include,  but  may  not  be  limited 
to,  all  7V-sohton  solutions.  In  the  process,  we  also  study  the  adjoint  equation  of 
the  eigenvalue  problem  for  the  BO  equation,  from  which,  for  real  potentials,  we  find 
a  very  simple  relation  between  the  fimctions  /?(A)  and  /(A),  introduced  in  Ref.[l]. 
Furthermore,  we  show  that  the  reflection  coefiicient  also  defines  a  phase  shift,  which 
can  be  interpreted  as  the  phase  shift  between  the  left  Jost  function  and  the  right 
Jost  function.  This  phase  shift  leads  to  an  analogy  of  Levinson’s  Theorem,  as  well 
as  a  condition  on  the  number  on  possible  bound  states  that  can  be  contained  in 
the  initial  data.  We  also  study  the  structure  of  the  scattering  data  and  the  Jost 
fimctions  for  pure  soliton  solutions,  and  obtain  remarkably  simple  solutions  for  these 
Jost  functions.  Since  they  are  examples  of  nongeneric  potentials,  they  demonstrate 
the  asymptotics  for  nongeneric  potentials.  We  then  carefully  detail  the  asymptotics 
in  the  limit  A  — »■  O’*"  for  generic  and  nongeneric  potentials.  Lastly,  we  show  how  to 
obtain  the  infinity  of  conserved  quantities  from  one  of  the  Jost  functions  of  the  BO 
equation,  and  also  how  to  obtain  these  conserved  quantities  in  terms  of  the  various 
moments  of  the  scattering  data. 
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I  INTRODUCTION 

The  Benjamin-Ono  (BO)  equation  [2] 


Ut  +  2uu,  +  Hu,,=0,  Hv{x)  =  -P  (1) 

where  P  indicates  the  principal  value  of  the  integral,  is  one  of  the  important  integrable 
equations  which  can  be  linearized  and  solved  by  an  inverse  scattering  transform  (1ST). 
It  has  important  physical  limits  of  internal  gravity  waves  in  a  stratified  fluid  [2]  as  well 
as  long  waves  in  a  stratified  shear  flow  [3,  4].  Thus  this  equation  is  very  important  for 
the  understanding  of  the  evolution  of  stratified  flows.  It  has  soliton  solutions  [5,  6,  7] 
and  a  Lax  Pair  [8,  9].  Interactions  of  BO  solitons  have  been  studied  by  Y.  Matsuno 
[10],  as  well  as  the  asymptotics  of  the  linearized  BO  equation,  which  describe  the  decay 
of  the  continuous  spectra.  A  perturbation  theory  for  the  1-soliton  solutions  of  the  BO 
equation  has  been  obtained  directly  by  Chen  and  Kaup  [11],  without  reference  to,  or 
the  use  of  the  1ST.  Stability  of  the  periodic  solutions  of  the  BO  equation  has  been 
studied  in  one  and  two  dimensions  [12] .  Physical  applications  of  soliton  solutions  of 
the  BO  equation  have  been  studied  in  various  contexts  by  Matsuno.  He  has  developed 
a  multisoliton  perturbation  theory  of  BO  solitons  [13],  and  has  applied  it  to  two  cases. 
These  cases  were  weak  dissipation  and/or  higher  order  dispersions,  and  he  determined 
the  shifts  in  the  soliton  amplitudes  and  the  phases  due  to  these  effects.  Effects  of 
uneven  bottoms  have  also  been  studied  by  Matsuno  [14]  wherein  he  has  demonstrated 
that  the  topography  could  capture,  repel,  or  even  generate  a  stationary  state. 

Fokas  and  Ablowitz  [1]  presented  the  1ST  solution  for  the  BO  equation.  This  solu¬ 
tion  was  novel  and  significant  in  several  ways.  First,  it  was  the  solution  of  the  direct 
and  inverse  scattering  problems  of  an  integro-differential  equation,  which  was  a  non¬ 
local  Riemann-Hilbert  problem  in  space.  As  such,  the  direct  problem  was  uniquely 
different,  requiring  Fredholm  theory  to  detail  the  properties  of  the  continuous  spec¬ 
tra,  the  discrete  spectra,  and  the  normalization  constants.  An  important  feature  of 
this  problem  is  that  the  generic  potentials  of  this  problem  will  in  general  only  vanish 
algebraic,  preventing  one  from  using  compact  support,  which  was  a  key  simplify¬ 
ing  feature,  useful  in  the  analysis  of  the  standard  Schrbdinger  and  Zakharov-Shabat 
eigenvalue  problems  [15],  Furthermore,  they  demonstrated  that  this  1ST  could  be 
formulated  without  having  to  detail  the  full  properties  of  the  scattering  data  of  the 
problem.  Also,  they  solved  these  problems  for  a  general  complex  generic  potential. 
The  direct  and  inverse  scattering  problems  for  a  related  scattering  problem  (origi¬ 
nating  from  a  different,  but  equivalent  Lax  pair),  has  been  studied  by  Coiffman  and 
Wickerhauser  [16].  There  they  have  described  many  of  the  features  of  this  scattering 
problem,  particularly  in  the  limit  of  very  small  initial  data. 

Here  we  extend  the  results  of  Ref.  [1]  in  several  ways.  First,  we  detail  the  symme¬ 
tries  of  the  scattering  data  when  the  potential  is  purely  real.  We  also  treat  the  case  of 
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nongeneric  potentials,  which  includes  all  the  iV-soliton  solutions.  As  we  shall  see,  the 
nongeneric  case  has  distinctly  different  asymptotics  from  the  generic  case  treated  in 
Ref.  [1].  When  the  potential  is  real,  one  can  also  relate  the  two  reflection  coefficients 
of  Ref.  [1],  by  a  complex  conjugation.  As  a  consequence  of  the  latter  result,  we 
then  find  that  the  inverse  scattering  equations  for  the  BO  equation  can  be  simplified. 
Then  by  including  the  asymptotics  of  the  Jost  functions  for  nongeneric  potentials, 
we  furthermore  find  that  the  inverse  scattering  equations  of  Ref.  [1]  can  be  extended 
to  include  this  case  as  well.  Also,  we  will  show  how  to  define  a  phase  shift  for  this 
scattering  problem,  the  change  of  which,  between  the  extreme  limits  of  the  spectral 
parameter,  can  be  related  to  the  number  of  bound  states.  This  phase  shift  is  related 
to  the  spectral  integral  in  the  Anderson-Taflin  conservation  law  [16,  17],  which  then 
shows  that  the  Anderson-Taflin  conservation  law  is  the  analog  of  Levinson’s  Theorem 
[18]  for  the  BO  equation. 

A  key  feature  of  these  new  results  is  the  treatment  of  the  adjoint  problem  of 
the  eigenvalue  problem  for  the  BO  equation,  by  which  we  are  able  to  define  inner 
products  of  the  Jost  functions.  Such  a  study  is  important  because  it  allows  one 
to  obtain  important  data  about  the  scattering  coefficients  and  the  orthogonality  of 
the  Jost  functions.  Such  a  treatment  will  undoubtedly  become  a  key  for  a  later 
development  of  a  general  perturbation  theory  for  the  BO  equation. 

In  the  next  section,  we  review  the  results  from  Ref.  [1].  They  are  reproduced  here 
only  for  completeness,  and  we  will  only  list  those  relations  that  we  will  need  to  refer 
back  to  and  use.  We  will  use  their  notation  for  the  definitions  of  the  Jost  functions  and 
the  various  scattering  coefficients.  In  later  sections,  we  will  introduce  new  quantities 
as  needed.  In  Section  III,  we  study  the  adjoint  problem  of  the  BO  eigenvalue  problem 
and  define  their  inner  products.  From  these  relations,  it  naturally  follows  that  the  two 
reflection  coefficients  of  Ref.  [1],  /(A)  and  f3{X),  are  related  by  a  complex  conjugation 
when  the  potential,  u{x)  is  purely  real.  Note  that  essentially  all  the  analysis  in  Ref. 
[1]  was  done  for  u{x)  complex  in  general,  in  which  case  these  two  reflection  coefficients 
would  be  independent.  Thus  it  is  not  surprising  that  if  one  requires  u{x)  to  be  purely 
real,  simplifications  will  occur.  In  Section  IV,  we  briefly  study  the  Jost  functions  for 
pure  AT-sohton  solutions.  These  jV-soUton  solutions  are  in  the  class  of  nongeneric 
potentials  and  have  features  distinctly  different  from  the  generic  potentials  studied 
in  Ref.  [1].  It  also  is  an  interesting  pedagogical  exercise  to  solve  the  BO  eigenvalue 
problem  for  these  pure  A^-soliton  solutions.  It  is  only  for  such  a  potential  that  closed 
form  solutions  for  the  Jost  functions  are  known.  We  shall  see  that  although  these  Jost 
functions  exist  for  almost  all  values  of  the  spectral  parameter.  A,  they  in  general  will 
have  logarithmic  singularities  in  the  upper-half  complex  x-plane.  These  singularities 
in  X  then  vanish  only  when  A  is  either  the  bound  state  eigenvalue  {Xj  <  0)  or  in  the 
continuous  spectra  (  A  real  and  >0).  In  Section  V,  we  extend  the  results  of  Ref. 
[1,  19]  for  the  asymptotics  of  the  Jost  functions,  in  the  limit  of  A  — >  O"*’,  to  include  the 
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general  case  of  nongeneric,  as  well  as  generic  potentials.  In  Section  VI,  we  simplify 
the  inverse  scattering  equations  of  Ref.  [1]  to  the  case  of  a  purely  real  potential, 
and  also  extend  them  to  include  the  possibility  of  nongeneric  potentials.  We  also 
specify  the  restricted  scattering  data  for  real  potentials  and  briefly  mention  the  time 
dependence  of  this  scattering  data  for  the  BO  equation,  which  is  given  in  Ref.  [1]. 
(For  simplicity  and  clarity,  we  shall  otherwise  ignore  the  time  dependence  in  all  the 
Jost  functions  and  the  scattering  data.)  In  Section  VII,  we  obtain  the  infinity  of  the 
conserved  quantities  for  the  BO  equation.  These  are  obtained  from  an  asymptotic 
expansion  of  the  Jost  function  W{x,X)  for  large  A.  We  also  give  these  conserved 
quantities  in  terms  of  the  scattering  data.  The  first  of  these  conserved  quantities 
relates  the  energy  in  the  BO  field  to  the  energy  in  the  scattering  data.  The  energy  of 
the  soUton  portion  and  the  energy  density  of  the  radiation  are  each  positive  definite 
and  additive.  From  the  conservation  of  the  total  area  of  the  BO  field  [16,  17],  we 
find  an  analogy  of  Levinson’s  Theorem  for  the  BO  equation.  Here  we  can  relate  the 
total  change  in  a  phase  shift  and  the  total  area  of  the  BO  field  to  the  total  number 
of  bound  states. 


II  BACKGROUND 


The  following  is  material  taken  from  Ref.  [1],  and  is  reproduced  here  for  completeness, 
since  we  shall  be  using  this  notation  and  will  need  to  refer  to  these  equations.  The 
eigenvalue  problem  for  the  BO  equation  may  be  taken  to  be 

Vx  —  iXv  —  ifitu]'''  =  0 ,  (2) 


where  v  is  the  eigenfunction,  analytic  in  the  upper  half  x-plane,  A  is  the  eigenvalue 
(spectral  parameter),  and  the  brackets  with  the  superscript “+”  indicates  that  one  is 
to  take  that  part  of  the  argument  which  is  analytic  in  the  upper  half  complex  x-plane. 
We  will  indicate  the  same  for  the  lower  half  ar-plane  with  a  superscript.  Thus,  as 
in  Ref.  [1],  we  take 

'  r 

^  2m 

for  e 


(3) 


^  —  x  =F  ze 

O'*".  There  are  two  Jost  functions  of  (2),  only  one  of  which  is  linearly  de¬ 
pendent.  First,  there  is  the  solution  N{x^  A)  defined  on  the  right,  and  there  is  the 
solution  M{x,X)  defined  on  the  left. 

N{x  ^  +00,  A)  ,  M(a:  ^  -oo,  A)  .  (4) 

These  solutions  exist  for  all  positive  values  of  A.  They  can  be  given  as  solutions  of 
integral  equations  [1],  such  as 

1  e~^^^dz 


1  6 

N(x,  1  I  ^ 


%e 


(5) 
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There  also  can  be  a  discrete  set  of  bound  state  solutions,  of  (2),  for  A  = 

Aj  <  0 ,  (j  =  1, 2, . . . ,  J)  where  J  is  the  total  number  of  such  bound  states,  possibly 
infinite.  Due  to  (2),  each  of  these  eigenstates  is  an  analytic  function  of  x  in  the  upper 
half  complex  x-plane.  The  bound  state  eigenfunctions  can  be  normalized  as 


+oo) 


1^ 

X 


(6) 


Upon  taking  the  limit  of  x  oo  of  (2),  from  (3)  and  (6),  it  follows  that 
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(7) 


and  one  finds  in  a  similar  manner  that  (6)  is  also  true  in  the  limit  oix  ^  — oo.  These 
functions  satisfy  the  homogeneous  form  of  (5),  which  is 

1  TOO  ^ 

I 

In  addition  to  these  eigenfunctions,  one  can  also  defined  two  auxiliary  functions, 
N{x,X)  and  M(x,  A),  each  of  which  are  a  particular  solution  of  the  inhomogeneous 
equation 

14  -  iXV  -  i[uV]-^  =  -iX ,  (9) 

and  are  defined  by  the  boundary  conditions 


N{x  +00,  A)  1 ,  M{x  -oo.  A)  1 .  (10) 

For  real  positive  A,  they  can  be  obtained  from  N  and  M  by  differentiating  with 
respect  to  A,  provided  /(A)  and  ^(A)  are  nonzero.  In  particular 

f{X)N{x,X)  =  e‘""a,[Ar(x,A)e-‘^1,  (11) 

g{X)M{x,X)  =  e'^"aA[M(x,A)e-‘^],  (12) 

where 

_ 1  roo 

_ 1  /*oo  _ 


They  also  satisfy  the  integral  equations, 

7f  (X.  A)e--  =  e-  -  I  £  A)  £  .  (15) 

These  functions  also  have  analytical  properties  with  respect  to  A.  M{x,X)  is 
analytic  in  the  upper  half  complex  A-plane  while  iV(x,  A)  is  analytic  in  the  lower  half 
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complex  A-plane.  By  Fredholm  theory,  Fokas  and  Ablowitz  then  were  able  to  show 
that,  in  the  appropriate  complex  A  half-plane,  as  A  approaches  any  of  the  bound  state 
eigenvalues,  each  of  these  functions  would  approach  the  limit: 

N{x,  A  ^  A,)  ^  +  {x  +  +  . . .  ,  (16) 


M{x,  A  \j)  ->  --!—^j{x)  +  {x  +  'rj)^j{x)  +  ....  (17) 

The  complex  quantity  jj  is  a  necessary  piece  of  the  scattering  data.  It  has  been  shown 
that  for  the  case  of  a  pure  7V-soliton  solution  [20]  and  also  indirectly  from  the  limit  of 
the  ILW  equation  [19],  its  imaginary  part  is  positive  and  is  given  by  —1/ (2Aj).  Later, 
in  Section  III,  we  shall  show  that  one  can  obtain  also  easily  obtain  this  same  result 
directly  from  the  inverse  scattering  equations  of  the  BO  equation. 

On  the  real,  positive  A-axis,  M{x,X)  is  related  to  N{x,X)  and  N{x,X)  by  the 
scattering  equation 

M{x,X)  =  N{x,X)+f3{X)N{x,X),  (18) 

where  P{X)  is  a  reflection  coefficient,  which  can  be  given  by 

/OO 

u{0M{^,X)e-^^d^.  (19) 

-OO 


Prom  the  above  scattering  equation  and  the  analytical  properties  of  N  and  M, 
Fokas  and  Ablowitz  obtained  a  key  relation  for  the  solution  of  the  inverse  scattering 
problem.  This  relation  relates  7i{x,X),  for  A  in  the  lower  half  complex  A-plane,  to 
the  eigenfunctions  N{x,  X)  on  the  real  A-axis  and  the  bound  states  $j(x).  This  is 


N{x,X)  =  1-1- 


p  p{X')N{x,X')dX' 
27ri  Jo  X'  —  X  +  ie 


Ai 


^jix) . 


(20) 


Requiring  (20)  to  satisfy  (9),  one  obtains  the  means  for  reconstructing  the  potential 
u{x)  from  the  scattering  data  and  the  Jost  frmctions.  This  is 

1  /•OO  _ _ _ 

=  e(X')N{x,X)dX'+i^ij{x),  (21) 


where  u{x)  =  [u|'*'(a:)-l-  complex  conjugate. 

The  first  of  the  inverse  scattering  equations  is  a  set  of  linear  dispersion  relations 
that  relate  the  bound  states,  ^j{x),  and  the  continuous  eigenstates,  N{x,X).  This 
relation  follows  from  (16)  and  (20),  upon  taking  the  limit  of  A  — ♦  Afc  in  (20).  Fokas 
and  Ablowitz  obtained 


{x  +  nfk)^k{x)  +iYj 

j^k 


p{X')N{x,  X')dX' 

A'-Afc 


=  1, 


(22) 
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which  is  Eq.  (27)  of  Ref.  [1]. 

In  addition  to  this  set  of  linear  dispersion  relations,  one  also  needs  a  linear  disper¬ 
sion  relation  for  the  continuous  spectra,  that  will  allow  one  to  also  determine  N{x,  A). 
This  will  follow  upon  integrating  (11)  with  respect  to  A,  and  then  using  (20)  to  elim¬ 
inate  N{x,X).  This  will  give  an  integral  equation  involving  integrals  with  respect  to 
A,  but  containing  only  N{x,X)  and  ^j{x),  which  is  Eqs.  (26)  and  (28)  in  Ref.  [1] 
and  Eqs.  (30)  -  (32)  in  Ref.  [19].  The  solution  of  these,  in  conjunction  with  (22), 
provide  the  solution  of  the  inverse  scattering  problem.  We  shall  not  give  these  forms 
here  since  they  were  obtained  for  complex,  generic  potentials  and  we  wish  to  simplify 
them  to  the  specific  case  of  real  potentials,  and  to  extend  them  so  that  they  may  also 
include  even  the  case  of  nongeneric  potentials. 


Ill  THE  ADJOINT  PROBLEM 

From  now  on,  we  shall  restrict  our  analysis  to  the  case  when  the  potential,  u{x),  is 
strictly  real.  This  will  provide  us  with  several  simplifications  and  will  relate  the  two 
Fokas  and  Ablowitz  reflection  coefficients.  Consider  the  adjoint  problem  of  (2).  It  is 

-b  iXv"^  —  iu[v'^]~  =  0 .  (23) 

Its  solution  is  rather  simple.  First,  take  the  part  of  (23)  which  is  analytic  in  the  lower 
half  complex  a:-plane.  Then  one  sees  that  one  simply  has  the  complex  conjugate 
of  (2).  Second,  take  the  part  of  (23)  which  is  analytic  in  the  upper  half  complex 
a;-plane.  This  is  a  simple,  first  order,  inhomogeneous  ODE.  Combining  these  two 
solutions  together,  one  has  that  a  solution  of  the  adjoint  problem,  N^{x,X),  is 

N\x,  A)  =  M\x,  X)-i  r  [uAr]+(e,  A)e-'^(*-«)de ,  (24) 

J  —  OO 

where  M  is  defined  by  (4).  The  last  term  in  (24)  is  the  part  which  is  analytic  in  the 
upper  half  complex  x-plane. 

We  define  the  adjoints  of  the  bound  state  eigenfunctions,  N^{x),  similarly. 

Nf{x)  =  3>;(a:)  -  i  T  .  (25) 

Next,  we  obtain  some  relations  between  our  Jost  functions  and  scattering  coef¬ 
ficients.  Note  that  N  and  M  are  not  independent  solutions.  Thus  they  must  be 
proportional  [19] .  For  A  >  0, 

M(a;,A)  =  r(A)N(x,A),  (26) 
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where  r(A)  has  two  equivalent  forms: 


r(A)  = 

roo 

l  +  i  w(e)M(e,A)e-‘^^de, 

•/-OO 

z= 

[l-</2'‘(0Af(?.A)e-‘«rfe]’‘. 

(27) 

Thus  by  (13)  and  (14), 

S(A)  =  r(A)/(A) . 

(28) 

Also,  from  (11),  (12),  (26)  and 

(28),  it  follows  that  for  real,  positive  A, 

(29) 

which  relates  F,  (3  and  g. 

We  will  need  the  asymptotics  of  these  various  Jost  functions.  These  all  follow 
from  the  definitions  and  (5),  (15),  and  (23).  For  a:  ^  +oo  and  A  >  0, 


while  for  x 


N{x,  A)  ->  , 

'N{x,X)  1, 

M{x,\)  e*^r(A), 

M{x,\)  l+y0(A)e‘^, 

N^{x,X)  e-*^^P(A), 

— oo  and  A  >  0, 


N{x,X) 

'N{x,X) 

M{x,  A) 
M{x,  A) 
N^ix,X) 


iXx 


r-^{x)e 

l-Me-, 


r(A) 


1, 

^-iXx 


For  the  bound  states,  for  x  — >  +oo  and  A  =  A_,-  <  0, 

X 

J—oo 

while  for  x  -*  -oo  and  X  =  Xj  <  0, 


i, 

X 

Nf{x)  ^  0. 


(30) 

(31) 

(32) 

(33) 

(34) 


(35) 

(36) 

(37) 

(38) 

(39) 


(40) 

(41) 


(42) 

(43) 
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Now,  from  the  Wronskian  relation,  we  can  relate  the  asymptotics  on  the  left  to 
those  on  the  right,  and  also  define  inner  products.  It  follows  that,  for  any  solution, 
v{x,X),  of  (2)  and  for  any  solution,  of  the  adjoint  problem,  (23), 


u^(a;,  X')v{x,  A) 


X) 

roo 

=  i{X  —  X')  /  v^{x,  X')v{x^  X)dx . 

J—oo 


(44) 


Thus  for  A'  =  X,v^  =  N'',v  =  N  and  from  (30),  (32),  (36)  and  (37),  it  follows  that 


r{A)r(A)  =  1 .  (45) 

Whence  we  may  define  a  real  phase  shift,  (9(A),  by 

r(A)  =  .  (46) 

Let  us  now  define  a  norm. 

A 

{v^\v)  =  j  v'^{x)v{x)dx .  (47) 

Then  from  the  Wronskian  relation,  (44),  it  follows  that 

(lV'‘(A')|Ar(A))  =  j^d(A-A'),  (48) 

where  5(A  —  A')  is  the  Dirac  delta  function.  Also 

(Ar/iyV(A))  =  0  =  (fV^(A')|$,> ,  (49) 

(Nj\it)  =0,  \(j  ^k,  (50) 

for  A  and  A^  real  and  positive. 

One  can  also  use  the  Wronskian  relation  with  the  function  iV  —  1 .  From  (9)  and 
(23),  one  obtains  the  general  relation 


7VV,A')[iV(a:,A)-l] 


la;=— oo 


TOO 

=  i(A-A')/  Af*(f,A')pV(f,A)-l|<i{ 

J  —<x> 

/oo 

^t/(e)[9v^(AO]-(Ode.  (51) 


We  consider  various  cases  for  (51).  For  A'  =  Xj  <  0,  from  (20),  (31),  (36),  (41),  (43), 

(49),  (50)  and  (51),  and  then  using  (7)  and  (25)  in  obtaining  the  second  line,  there 
results 


(JV/14.,)  =  i  r  u{i)[N}\-(i)d^, 

J—OO 

=  — 27rAj . 


(62) 
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This  gives  us  our  nonzero  inner  products  of  the  bound  states.  Next,  for  X'  >  0,  from 
(20),  (31),  (34),  (36),  (39),  (48),  and  (51),  and  again,  then  using  (13),  (24),  (26)  and 
(45)  in  obtaining  the  second  line,  there  results  the  key  relation,  necessary  for  relating 
/(A)  and  /?(A). 


/3(x')  =  -ir{x')  r 

J  —  oo 

=  -2iriX'r{X'), 


from  which  we  have 


m  = 


(53) 


(54) 


27riA  ’ 

which  is  exactly  the  relation  that  one  would  find  from  the  linear  limit.  When  we 
combine  (29)  with  (28)  and  (53),  we  also  find 


dx9{X)  = 


iirX  ' 


(55) 


Thus  the  phase  shift  will  always  be  a  monotonic  function  of  A. 

We  can  now  demonstrate  that  the  imaginary  part  of  7^-  is  related  to  the  eigenvalue 
Aj-,  as  mentioned  under  Equation  (17).  As  mentioned  before,  this  already  has  been 
shown  to  follow  from  the  BO  limit  of  the  ILW  equation  [19].  Here,  we  simply  show 
how  to  obtain  this  result  purely  from  the  BO  eigenvalue  problem.  Take  the  inner 
product  of  (22)  with  the  state  N^{x),  and  then  take  the  imaginary  part.  One  obtains 


Im  f  Nf{x)[x^j{x)  —  \\dx  =  27rAjIm7j 

J  —  00 


(56) 


Now  use  (25).  The  second  part  of  N^[x)  is  the  part  of  N^(x)  which  is  analytic  in 
the  upper  half  complex  x-plane,  as  is  also  [a;4>j(a;)  -  1],  with  the  latter  vanishing  like 
l/x  for  large  |a;|.  Thus  this  part  gives  a  zero  contribution  to  the  integral.  The  first 
part  of  Nf{x)  combined  with  the  x^j{x)  term  gives  a  real  quantity,  so  neither  does 
it  contribute.  What  one  is  left  with  is  simply 

/oo 

^j{x)dx  =  27rAjIm7,- .  (57) 

-OO 

Evaluating  the  final  integral  along  an  infinite  semicircle  in  the  lower  half  complex 
x-plane,  gives  the  final  result 

Im7i  =  ^  •  (58) 
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IV  JOST  FUNCTIONS  FOR  A  PURE  SOLITON 

SOLUTION 

Although  AT-sohton  solutions  of  the  BO  equation  have  been  studied  often,  so  far  it 
seems  that  there  has  been  no  study  of  the  dost  functions  for  a  single  soliton  Thlse 
are  very  simple,  very  instructive  solutions  and  their  structure  demonstrates  some  of 
the  features  of  this  scattering  problem.  Thus  we  give  them  here  and  point  out  some 
o  their  more  important  and  unique  features. 

When  the  potential  is  a  pure  1-soliton  solution,  the  reflection  coefficient  /?(A)  must 

ha^(A)  =  0  by  (54).  We  can  immediately  obtain  the  bound  states 
from  (22)  and  then  N{X)  from  (20).  Let  the  bound  state  eigenvalue  be  A,  =  -la 
where  o  >  0.  Prom  (22),  it  follows  that  ^ 


x  +  'ri 


(59) 


where  from  (58),  we  take 

7i  =  +i/a, 

where  Xo  will  be  the  center  of  the  soliton  and  a  will  be  the  amplitude.  Note  that  since 

n  >  0,  we  have  $1(1)  analytic  in  the  upper  half  complex  i-plane.  The  potential 
now  be  recovered  by  (21),  which  gives 


can 


uY'{x)  = 


a:  +  7i  ’ 


u{x)  = 


2/a 


(61) 


(i  +  7i)(H-7;)  ' 

Although  we  are  able  to  recover  the  iV-soliton  potential,  at  the  moment,  we  have 
no  mecham^  for  the  direct  construction  of  the  Jost  function,  N{x,X).  Clearly,  we 
can  obtain  N{x,X)  from  (20)  since  ^(A)  =  0.  However,  since  /(A)  is  zero,  we  can 
only  obtain  from  (11),  the  seemingly  trivial  result 


N{x,X)  =  n{x)e 


iXx 


(62) 


where  n(a:)  is  independent  of  A.  This  result  will  hold  for  any  of  the  Af-soliton  solutions 
As  we  shall  see  later,  one  can  determine  n{x),  but  it  will  be  necessary  to  first  establish 
tne  asymptotics  of  the  nongeneric  case  in  the  limit  of  A  0"^.) 

Otherwise,  the  only  means  that  we  have  at  the  moment  for  reconstruct  of  Nix  A) 
IS  by  means  of  the  direct  scattering  problem,  (2)  and  (4) .  We  illustrate  how  to  do  this 
in  general  for  the  one-soliton  potential,  (61).  The  general  multisoliton  case  can  be 
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treated  in  the  same  way.  Noting  that  N  is  to  be  analytic  in  the  upper  half  complex 
x-plane,  the  function  can  be  evaluated  using  the  residue  theorem,  which  gives 


[uTV]'*'  =  u{x)N{x,  A)  +  i 


N{xo  +  i/a,  A) 

X  —  Xo  —  i/o, 


(63) 


Note  that  the  last  term  has  a  pole  in  the  upper  half  complex  x-plane.  It  now  is 
relatively  easily  to  integrate  (2)  under  the  boundary  condition  (4).  The  result  can  be 
expressed  in  the  form 


N{x,  A)  =  e’ 


iXx 


X 


|[l+W{i„+j,A){(l+^)eV“S,(z)  ■ 


X-XO  + 


a  X 


Xq 


— }] .  (64) 


where  z  =  i\(x  —  Xq  —  if  a)  and  Ei{z)  is  the  exponential  integral,  defined  by 


fOO 

Ei{z)^I^  —dt,  larg(z)|<7r.  (65) 

The  expressions  (63)  and  (64)  still  contain  the  unknown  factor  N{xo  +  i/a,X).  How¬ 
ever,  in  view  of  the  fact  that  Ei{z)  has  a  logarithmic  singularity  at  z  =  0  (  jEi(2r)  ~ 
—'y—\Q.z-\-0{z)  as  z  — >■  0,  where  7  is  Euler’s  constant),  which  is  in  the  upper  half  com¬ 
plex  x-plane,  one  can  see  that  in  general,  the  result  (64),  is  not  analytic  in  the  upper 
half  plane,  unless  the  products  of  the  coefficients  of  this  term,  N{xQ^-i/a^  A)(l  -t-2A/a) 
happen  to  vanish. 

E  we  choose  N{xo  +  i/a,  A)  =  0,  then  (64)  takes  on  the  particularly  simple  form 


N{x,  A) 


^i^x^-xo^-i/a 
X  —  Xo  +  i/a  ’ 


(66) 


which  is  exactly  of  the  form  predicted  by  (62).  Furthermore,  in  the  limit  of  a:  — —00, 
from  (4),  (26),  and  (66),  it  follows  that  r(A)  —  1  for  the  soliton  solution.  Note  that 
this  is  in  accordance  with  (45),  (46)  and  (55),  since  9  must  be  a  constant  for  /?(A)  =  0. 
It  also  can  be  confirmed  directly  from  (27),  (61)  and  (66). 

The  solution  (64)  is  the  continuous  spectrum,  A  >  0.  For  the  bound  states,  we 
must  have  A  <  0.  This  solutions  will  be  the  same  as  (64),  except  that  the  first  term, 
1,  in  the  first  bracket  of  (64)  will  be  absent.  Again,  the  term  Ei{z)  can  be  made  to 
vanish,  but  now,  by  taking  (1  -f  ^)  =  0.  However,  this  is  just  a  restatement  of  (58), 
in  that  the  imaginary  part  of  71,  (1/a),  is  related  to  the  bound  state  eigenvalue  by 
Ai  =  —\a.  With  this,  the  solution  reduces  to  the  bound  state  solution  given  by  (59). 

Another  important  point  of  these  nongeneric  solutions  is  that  from  (61)  and  (66), 
one  obtains 

/OO 

u{x)N{x,\)e-^^^dx  =  Q.  (67) 

*00 

As  we  shall  shortly  see,  this  has  implications  for  the  asymptotics  of  the  Jost  functions 
for  nongeneric  potentials  in  the  limit  of  A  — >  O'*’.  It  shows  that  the  nongeneric 
potentials  can  have  this  integral  as  being  zero,  in  contrast  to  generic  potentials,  for 
which  this  integral  is  always  nonzero  [1,  19]. 
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V  ASYMPTOTICS  OF  THE  JOST  FUNCTIONS 


In  order  to  extend  the  inverse  scattering  equations  of  Ref.  [1]  to  include  the  nongeneric 
cases,  it  will  be  necessary  to  address  the  asymptotics  of  these  Jost  functions,  N{x,  A) 
and  N{x,  A),  and  the  reflection  coeflficient,  /3{X),  in  the  limit  of  A  — ^  O"''.  Let  us  start 
with  the  Green’s  function  form  of  the  solutions  [1],  (5)  and  (15).  One  can  reduce 
these  to: 

fOO 

N{x,  X)e-^^  -l+i  u{ON{^,  A)e-‘^^d^ 

Jx 

1 

N{x,  -e-‘^  +  i  r  u(()N{i, 

Jx 

1  fOO  _ 

-^J_^<ON{^,X)e-^^^Ei[iX{x-^-i€M  =  0, 

where  Ei{z)  is  the  exponential  integral,  defined  by  (65).  Now  define 

^i{z,  A)  Ei^Xi^z  —  ie)]  +  7  "b  (^'^)  >  (”^6) 

where  7  is  Euler’s  constant.  In  the  limit  of  A  0"^,  we  have 

£i(z,A^0+)  ^ -lnz  +  0(zA).  (71) 

The  asymptotics  given  in  [1]  and  [19]  handle  the  generic  case  quite  well.  The 
asymptotics  for  the  nongeneric  case  can  also  be  included  in  a  like  manner,  provided 
one  allows  for  the  vanishing  of  an  integral,  namely  (67),  in  the  limit  of  A  0+.  To 
do  this,  first  note  that  without  any  approximations,  (68)  can  very  nicely  be  put  into 
the  form 

/*oo 

Nq{x,  A)e-‘^  -l+i  u{^)No{^,  A)e-‘^^d^ 

Jx 

1 

where 

«o(i,A)  =  N(x,\)IU{\), 

H(\)  =  \-[l  +  \n{i\)\+j"^u(ON{i,\)e-‘^di, 

_ _ 1 _ 

1  +  [7  +  In  (iA)|i  uiONoi^,  X)e-^^d^  ’ 

with  J\f(X)  being  a  constant,  independent  of  x.  As  is  clear  from  (72),  the  function 
No{x,X)  satisfies  a  nonhomogeneous  integral  equation  and  therefore  is  uniquely  de¬ 
termined.  Thus  it  cannot  be  normalized.  Furthermore,  for  A  0"^,  due  to  (70),  the 


(72) 

(73) 

(74) 


(68) 

(69) 
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solution  for  No{x,X)  can  be  expanded  in  an  asymptotic  power  series  about  A  =  0, 
and  when  so  expanded,  will  have  no  In  A  terms  present.  Thus 


Nq{x,  a)  — ^  Noo{x)  +  0(A) ,  (75) 

where  Noo{x)  is  the  solution  of  (72)  for  A  =  0  and  satisfies 

roo  1  roo 

Noo{x)-l+iJ^  u{ONoo{Od^  +  ^J  ^u{ONoo{O^ri{x-^-ie)d^  =  0.  (76) 

Clearly,  N{x,  A)  will  have  a  similar  structure  and  can  be  put  into  a  similar  form. 

_  1  _  p—iXx  i*oo  _ 

N„{x,  A)e-“*  -  1  +  +  i  u(f)W„K, 

-^/_”u(fl^c.teA)e-“f<J?£,[iA(3;-{-ie)|  =  0,  (77) 

where 


A^o(a:,A) 

77(A) 


iV(a;,A)/:A7(A), 

1 


1  +  [7  +  ln(iA)]i/_~  n(OiVo(^,A)e->^de  ' 


(78) 

(79) 


However,  these  asymptotics  will  be  slightly  different.  Due  to  the  presence  of  Ar(A)  in 
(77),  we  have  instead 


W„(i.A)-.W„„(x)+0(A)+0(^),  (80) 

where  the  function  Noo{x)  is  the  same  fimction  as  in  (75),  which  is  given  by  (76). 

The  terms  A7(A)  and  77(A)  contain  the  asymptotic  forms  of  the  A  dependence  of 
the  dost  functions  for  all  possible  cases,  generic  and  nongeneric.  We  already  know 
that  there  are  solutions  for  which  the  integral 

/oo=  r  uiONooim.  (81) 

J—OO 

is  exactly  zero,  (67),  which  are  the  AT-soliton  solutions.  It  is  not  known  if  there 
are  other  examples,  but  it  is  possible  that  initial  value  problems  where  u{x)  has  a 
vanishing  small  amplitude  with  a  vanishing  area  may  also  have  this  integral  zero. 

Before  continuing,  we  briefly  comment  concerning  the  asymptotic  forms  of  these 
dost  functions.  It  should  be  clear  that  the  above  forms  will  be  useful  only  as  long  as 
|Ax|  <<  1,  since  it  is  only  then  that  we  can  And  the  expansion  (71)  useful.  However 
all  equations  given  here,  from  (68)  to  (80),  are  exact,  except  for  the  three  asymptotic 
forms;  (71),  (75)  and  (80).  Now,  take  a  wider  view  of  these  asymptotics.  For  any 
given,  sufficiently  small  value  of  A  >  0,  there  will  always  be  4  diff’erent  regions,  each 
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containing  different  asymptotics  for  the  Jost  functions.  First,  there  is  the  region 
where  u(x)  is  essentially  nonzero.  Second,  (and  considering  only  the  regions  to  the 
right  of  the  potential,  there  are  also  similar  regions  to  the  left)  there  is  the  region 
where  u(x)  has  become  effectively  zero,  but  one  still  has  x\  «  1.  Here,  from  (76) 
one  has  Nm{x)  ->  1  -  (In  x)^/oo,  which,  when  the  integral  is  nonzero,  contains 
the  Inx  dependence  found  in  |1|  and  119|  for  the  generic  case.  Continuing,  for  any 
given  A  >  0  there  wiU  exist  values  of  x  for  which  Ai  =  0(1),  with  the  potential 
also  being  effectively  zero.  This  is  the  third  region,  where  we  no  longer  can  use  the 
asymptotics  provided  by  (70),  but  must  return  to  the  original  form  given  by  (68) 

IS  IS  a  trmsition  region,  where  the  nature  of  the  solution  is  changing  from  an 
mterior  solution  to  an  exterior  solution.  The  asymptotics  here  are  now  determined 
by  B.(xA)  instead.  Lastly,  there  is  the  fourth  region,  where  for  any  given  A  >  0  rr 
has  become  so  large  that  xA  »  1.  Here  there  is  no  remaining  Inx  dependence  in  the 
dost  Actions  and  instead  the  Jost  functions  are  approaching  the  appropriate  values 

jO  —  OO. 

Now,  let  US  turn  to  the  asymptotics  of  the  scattering  data.  As  in  Ref.  [1|  we 

can  obtain  tbs  information  from  the  spatial  integrals  over  a  product  of  the  potential 
with  a  Jost  faction.  To  evaluate  the  asymptotics  of  such  a  integral,  we  only  ne ^ 
o  consider  the  i^ymptoti.  of  the  Jost  function  in  the  first  region,  smce  we " 

bZ"rvah7S " : 

IT «8  Je  ^  ■  tr  1  r  ('8),  and  (27).  We  then 

relations  ^  ™s  gives  us  the  exact 


m 

1 

r(A) 


/OO 

w(0A^(e,A)rfe, 

■OO 

/OO  _ 

/OO 

^«(f)JV(^,A)e'‘«d^. 


(82) 

(83) 

(84) 

to3  yl)  fi  d  V  ‘he  asymp. 

nonzero.  If  it  is  nonzero,  then  we  have  the  generic  case  where 


Ar(A) 
A7(A) 
JV  (x,  A) 


27r  1  V 


27r 


1 


/ooln  A 
2-kNoo{x) 
Iqq  In  A 


(85) 

(86) 

(87) 
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N{x,  A) 

m 

m 


2TrNoo{x) 


Iqo  In  A 

27r  1 

+  0( 


In^A' 


In  A 


In^A 


), 


)■ 


(88) 

(89) 

(90) 


If  the  integral  Iqq  is  zero,  we  then  have  the  nongeneric  case,  where  instead 


J\f(^X)  l+0(AlnA),  (91) 

77(A)  — ^  l  +  0(AlnA),  (92) 

N{x,X)  Noo{x)  +  0{XlnX),  (93) 

TV  (a:,  A)  TVoo(^)  +  0(Aln  A) ,  (94) 

/?(A)  O(AlnA),  (95) 

r(A)  ^  l+0{AlnA).  (96) 


Prom  (55),  (90)  and  (96),  we  now  have  the  full  solution  for  ^(A).  It  is 


Given  initial  scattering  data  for  the  BO  equation,  it  is  a  simple  matter  to  determine 
whether  one  has  the  case  Iqq  zero  (the  nongeneric  case)  or  nonzero  (the  generic  case). 
Simply  look  at  the  asymptotics  of  ^(A  0+)  and  determine  how  fast  the  reflection 

coefficient  is  vanishing.  Note  that  we  do  not  assume  that  the  nongeneric  case  only 
includes  the  N-soliton  solutions,  which  requires  p{X)  —  0.  Although  the  nongeneric 
case  does  contain  the  A-soliton  solutions,  we  have  not  been  able  to  exclude  the 
possibility  of  the  existence  of  nongeneric  cases  where  /?(A)  ^  0. 


VI  SIMPLIFICATION  OF,  AND  EXTENSION  OF 
THE  1ST  FOR  THE  BO  EQUATION 

With  the  above,  we  may  now  simplify  the  last  linear  dispersion  relation  (Eq.(26)  of 
Ref.  [1])  to  the  case  of  a  real  potential,  and  also  extend  it  so  as  to  include  nongeneric, 
as  well  as  generic  potentials.  We  shall  proceed  as  in  Ref.  [1]  except  for  one  key 
difference.  We  shall  first  subtract  off  a  collection  of  terms  that,  upon  integration, 
would  give  singular  integrals. 

Consider  (11)  and  (20).  If  one  attempts  to  integrate  the  resulting  expression 
directly,  one  obtains  a  singular  result  in  the  limit  of  A  0+.  In  Ref.  [1],  this  was 
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handled  by  subtracting  off  the  singular  part  of  the  reflection  coefficient,  /(A),  into  the 
factor  /s(A),  and  then  integrating  the  terms  in  such  a  fashion,  so  that  the  singular 
integrals  would  can  cell.  Here,  we  shall  use  a  slightly  different  procedure  for  handhng 
these  singular  terms.  We  shall  simply  subtract  off  the  limit  of  a  Jost  function  in  the 
limit  of  A  — >  O'*'.  This  Jost  fimction  vanishes  in  this  limit  for  generic  potentials,  but 
for  nongeneric  potentials,  it  is  a  key  component  of  the  solution.  To  see  this,  take  the 
limit  of  A  0+  in  (20).  There  results 


,  1  mX')N(x,X')d\' 


(98) 


As  we  have  seen  in  the  previous  section,  the  function  A^(a;,  0)  is  either  0  or  Noo{x), 
depending  on  whether  Iqq  is  nonzero  (generic)  or  zero  (nongeneric).  Note  that  in 
either  case,  the  integral  in  (98)  is  not  divergent  as  A  — ^  O'*".  (Note  also  that  in  the 
generic  case,  (98)  is  then  an  identity  and  therefore  a  constraint  on  the  Jost  functions 
and  scattering  data,  particularly  the  parts  that  give  rise  to  the  singular  integrals. 
Thus  this  proceedure  could  have  been  used  equally  well  in  Ref.  [1],  which  would  have 
given  them  a  simplify  version  of  their  Eq.  (26).) 

Let  us  now  use  (98)  to  subtract  off  the  parts  of  (20)  which  would  give  singular 
integrals.  Simply  subtracting  gives 


N{x,X)=N{x,0)  + 


p{X')N{x,  X')dX' 
2'Ki  Jo  A'(A'  —  A  +  ie) 


1 

A,(A-A,) 


(99) 


Note  that  the  last  two  terms  are  now  proportional  to  A. 

Using  (11)  to  eliminate  N{x^X)  from  (99),  then  as  in  Ref.  [1],  we  may  integrate 
the  result  from  A  ^  O"*"  to  A,  and  noting  that  N{x,0)  =  N{x,0),  we  obtain 


_  _  1 
N{x,  A)e-‘^*  =  N{x,  0)  -  iw{x,  A,  0)N{x,  0)  -  ^  — $j(x)w(x,  A,  Xj) 

i=i  ^3 

^P(X')N{x,  X')w(x,  A,  A') ,  (100) 


27r 


A' 


where 


»(x,A.A)  =  -y^ 


(101) 


for  X'  real  (and  either  positive  or  negative).  This  is  the  simplified  and  extended  form 
of  Eqs.  (26)  and  (28)  of  Ref.  [1].  As  one  may  readily  verify  from  (87)-(95),  all 
the  integrals  in  (100)  and  (101)  are  well  defined  and  nonsingular  at  the  lower  limit, 
except  for  the  term  iy(a;,A,  0)  in  the  generic  case  {Iqq  ^  0).  However  in  this  case, 
since  N{x,0)  is  zero  anyway,  the  quantity  u;(x,  A,0)  is  never  required. 

In  the  nongeneric  case,  (/qo  =  0),  we  close  our  set  of  equations  by  including  (98)  in 
the  set  of  linear  dispersion  relations.  Our  final  set  of  inverse  scattering  equations,  valid 
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only  for  real  potentials,  but  extended  to  include  nongeneric  potentials,  are  therefore 
Eqs.  (22),  (98),  and  (100),  with  recovery  of  the  potential  being  achieved  by  (21). 
Note  that  once  these  solutions  are  obtained,  then  one  may  construct  N{x,X)  from 
either  (20)  or  (99). 

For  iV-soliton  solutions,  since  /?(A)  =  0,  all  terms  in  (100)  vanish,  except  for  the 
first  term,  N{x,0),  which  can  now  be  given  by  (98).  In  this  case,  (100)  reduces  to 
(62)  with  its  value  being  given  by  (98). 

This  particular  structure  of  a  set  of  inverse  scattering  equations  deserves  some 
comments.  As  already  noted,  in  the  generic  case  (/qo  ^  0),  (98)  gives  a  constraint 
on  the  scattering  data,  since  the  left-hand  side  is  zero.  This  can  be  understood  as 
follows.  Due  to  (87),  (88),  and  (89),  the  reflection  coefficient  and  the  Jost  functions 
have  been  predetermined  at  one  point,  A  — »  O'*".  Thus  one  degree  of  freedom  is  missing 
from  the  continuous  spectrum,  which  is  taken  up  by  the  nonzero  value  of  Iqq.  On  the 
other  hand,  in  the  nongeneric  case,  (/qo  =  0),  then  the  reflection  coefficient  and  the 
Jost  functions  are  not  predetermined  in  this  limit  (only  the  order  is  known,  not  the 
coefficient),  and  (98)  and  (100)  then  are  independent  equations,  with  no  degrees  of 
freedom  missing. 

The  time  dependence  of  the  scattering  data  for  the  BO  equation  was  first  given 
in  Ref.  [1].  We  include  it  here  for  completeness,  and  it  is  simply 


dXj{t) 

dt 

dljjt) 

dt 

dt 


0,  J  =  l,2 . J, 


(102) 

(103) 

(104) 


VII  AN  INFINITY  OF  CONSERVATION  LAWS 

Although  it  has  long  been  known  that  the  BO  equation  has  an  infinity  of  conserved 
quantities  [8,  9,  21],  no  general  scheme  has  ever  been  given  for  their  generation  from 
the  scattering  problem  [22,  23,  24],  nor  has  anyone  ever  given  them  in  terms  of  the 
scattering  data.  This  we  shall  do  here. 

Start  with  the  eigenvalue  equation  for  77{x,X),  (9),  and  its  time  evolution 

-  iX'N  -  i[u77]+  =  -iX ,  (105) 

Nt  -  2XN^  -  -  2[r]+7V  =  0 .  (106) 

Then  from  (1),  (105)  and  (106),  one  has  that 


[u{x)N(x,  A)|fda;  =  0  . 


(107) 
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Thus  the  quantity  /_^  u{x)N{x,  X)dx  is  conserved.  This  conserved  quantity  can  be 
evaluated  in  terms  of  the  scattering  data.  It  follows  from  (20)  that 


j=l  ^ 


(108) 


(109) 


where  (7),  (13),  (M)  and  (104)  have  been  used. 

Now,  expand  N(x,X)  in  inverse  powers  of  A,  as 

3v(x.A)  =  f;(zl)::^i±iM. 

Substitute  this  into  (105).  One  can  now  determine  N^+iix)  successively  from; 

■^1  —  1  ,  Nn+i  =  iNn,x  +  [wA^n]"^,  (n  >  1)  .  (110) 

Expanding  (108)  also  in  an  asymptotic  series  for  A  large,  then  substituting  (109)  into 
the  result  and  comparing  the  coefficients  of  A'”  on  both  sides,  we  obtain 

/oo  _ 

u{x)Nn{x)dx , 

oo 

1  9  J 

=  nmy)dy  +  2^j:(-Xjr-',  (n>2).  (Ill) 

i=i 

winch  are  the  inMty  of  conserved  quantities.  Qearly,  one  can  generate  ]V.(x)  to  any 
order  and  therefore  any  conserved  quantity.  FWthermore,  note  that  each  conserved 
1  y  IS  smp  y  the  (n  2)th  moment  of  the  radiation  density,  0'{\)/3(X)  plus 
a  corsespondrng  term  for  the  soliton  contribution.  We  note  that  L  fiLtlL 

quanutyj  i2>  is 


/oo  2  j  - 

^  =  2wg(-A,.)  +  _  £° dX0-{X)0(X) , 


(112) 


wbch  relate  the  energy  in  u(x)  to  the  energy  in  the  scattering  data. 

from  the'  b  '  conserved,  its  expansion  cannot  be  obtained 

from  the  above  Instead,  .t  has  to  be  obtained  by  a  separate  consideration,  which  we 

rt!  tt  h  f  r  (98)  by  n(x)  and  integrating 

irom  oo  to  oo,  by  which  we  obtain 

O  T  1  ^  roo 

^  ~2TrL  ('')/^(-')  +  /^“(i)dx,  (113) 

which  relates  the  number  of  bound  states,  J,  and  the  total  “quanta”  of  radiation  to 
the  total  area  under  u(x).  This  is  the  Anderson-Taflin  |17J  conservation  law  and  its 
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An  important  feature  of  this  conservation  law  follows  from  (97).  Expressed  it  as 

/OO 

u{x)dx ,  (114) 

■OO 

shows  that  the  total  number  of  bound  states  depends  only  on  the  value  of  the  to¬ 
tal  phase  shift,  ^(oo),  and  the  total  area  under  the  profile,  both  of  which  are  time 
independent.  Also,  we  see  that  for  a  given  area,  there  is  no  limit  on  the  number  of 
solitons  that  could  be  generated,  unless  the  total  phase  shift  was  somehow  limited. 
But  on  the  other  hand,  from  (112),  we  see  that  for  a  given  total  area  of  a  profile  and 
a  given  total  energy,  there  is  a  limit  on  the  total  number  of  possible  solitons  that  are 
larger  than  any  given  size. 

The  form  of  (114)  is  very  suggestive.  Since  both  J  and  the  total  phase  shift  are 
positive  definite,  we  see  that  the  total  area  under  the  profile  is  expressible  as  a  positive 
quantity,  27rJ,  minus  another  positive  quantity,  0(oo).  Thus,  in  some  fashion,  it  seems 
that  we  can  expect  the  positive  regions  of  u{x)  to  be  the  sources  of  the  solitons  while 
the  negative  regions  would  be  the  sources  of  the  radiation.  In  fact,  in  the  small 
dispersion  limit  of  the  BO  equation,  one  can  derive  exactly  this  condition  [25,  26,  27], 
whereby  the  strictly  positive  parts  of  the  profile  area  determine  the  number  of  solitons. 


VIII  CONCLUSIONS 


We  have  restricted  the  results  of  Ref.  [1]  to  real  potentials,  obtaining  a  reduction 
in  the  number  of  reflection  coefficients  required.  We  have  also  extended  the  same  to 
include  nongeneric  potentials  and  have  generalized  their  inverse  scattering  equations 
to  include  this  case. 

When  the  potential  is  real,  we  have  shown  that  /(A)  is  directly  proportional 
to  P*{X).  We  have  also  generalized  the  asymptotics  of  the  dost  functions  and  the 
scattering  data  to  the  nongeneric  case.  With  only  one  minor  change,  we  have  modified 
the  set  of  inverse  scattering  equations  introduced  in  Ref.  [1],  by  the  introduction  of 
one  new  function,  N(x,0),  which  vanishes  for  generic  potentials,  but  is  nonzero,  and 
is  a  key  function,  for  nongeneric  potentials. 

We  note  that  although  A^-soliton  solutions,  for  which  /?(A)  vanishes  identically,  are 
one  set  of  examples  of  nongeneric  potentials,  we  do  not  know  if  there  are  other  sets. 
This  class  of  potentials  is  distinguished  from  generic  potentials  by  the  vanishing  of  a 
single  integral,  /qq,  complex  in  general,  and  thus  with  only  two  degrees  of  freedom. 
Considering  this,  it  then  seems  quite  unlikely  that  the  vanishing  of  this  single  integral, 
with  only  its  two  degrees  of  freedom,  would  be  sufficient  to  guarantee  that  /?(A),  would 
have  to  vanish  identically,  for  all  A.  Considering  such,  then  one  expects  that  there 
must  exist  other  examples  of  nongeneric  potentials  besides  A^-soliton  solutions.  We 
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have  already  noted  that  one  such  example  could  be  a  zero  area  pulse  with  a  vanishing 
small  amplitude. 

We  have  also  presented  the  necessary  recurrence  relations  for  generating  the  infin¬ 
ity  of  conserved  quantities  from  the  Jost  function,  N{x,  A),  and  also  have  given  these 
quantities  in  terms  of  the  scattering  data.  We  have  shown  that  the  Anderson-Taflin 
[17]  conservation  law  is  the  analogy  of  Levinson’s  Theorem  for  the  BO  equation  and 
that  it  relates  the  total  area  of  the  BO  field  to  a  total  phase  shift  and  the  total  number 
of  bound  states. 

We  have  also  obtained  the  inner  products  of  the  Jost  functions  with  their  adjoints. 
These  functions  can  be  expected  to  be  a  closed  set  for  any  function  analytic  in  the 
upper  half  x-plane,  and  a  proof  of  such  could  be  expected  to  follow  from  a  similar 
proof  for  the  ILW  equation  in  Ref.  [28].  If  we  adjoin  to  these  functions  their  complex 
conjugates,  then  the  total  set  would  be  expected  to  be  complete  for  any  L2  function. 
Also,  from  the  inner  products,  assuming  closure,  one  finds  a  very  suggestive  structure 
of  the  closure  relation.  It  may  be  possible  to  interpret  this  in  some  manner,  perhaps 
equivalent  to  bound  state  eigenvalues  being  given  by  the  zeros  of  r(A),  somehow 
continued  into  the  negative  real  A-axis.  This  was  a  key  feature  of  the  Zakharov-Shabat 
1ST  [15].  This  also  seems  to  be  implied  from  the  BO  limit  of  the  ILW  equation  [19]. 


ACKNOWLEDGMENTS 

The  authors  gratefully  acknowledge  valuable  comments  by  A.S.  Fokas,  and  also 
those  of  an  anonymous  referee.  This  research  was  supported  in  part  by  the  Office  of 
Naval  Research,  under  Grant  No.  N00014-95-1-0323. 

Acknowledgement  and  Disclaimer:  Effort  sponsored  in  part  by  the  Air  Force  Office 
of  Scientific  Research,  Air  Force  Materials  Command,  USAF,  under  contract  number 
F49620-96-C-0031 .  The  US  Government  is  authorized  to  reproduce  and  distribute 
reprints  for  governmental  purposes  notwithstanding  any  copyright  notation  thereon. 
The  views  and  conclusions  contained  herein  are  those  of  the  authors  and  should 
not  be  interpreted  as  necessarily  representing  the  official  policies  or  endorsements, 
either  expressed  or  imphed,  of  the  Air  Force  Office  of  Scientific  Research  or  the  US 
Government. 


References 

[1]  A.  S.  Fokas  and  M.  J.  Ablowitz,  Stud.  Appl.  Math.  68,  1  (1983). 


21 


[2]  T.  B.  Benjamin,  J.  Fluid  Mech.  29,  559  (1967);  R.E.  Davies  k  A.  Acrivos,  J. 
Fluid  Mech.  29,  593  (1967);  H.  Ono,  J.  Phys.  Soc.  Japan  39,  1082  (1975). 

[3]  S.  A.  Maslowe  and  L.  G.  Redekopp,  J.  Fluid  Mech.  101,  321  (1980). 

[4]  R.  Grimshaw,  Stud.  Appl.  Math.  65,  159  (1981). 

[5]  K.  M.  Case,  Proc.  Nat.  Acad.  ScL,  USA  75,  3563  (1978);  76,  1  (1979). 

[6]  H.  H.  Chen,  Y.  C.  Lee,  and  N.  R.  Pereira,  Phys.  Fluids  22,  187  (1979). 

[7]  Y.  Matsuno,  J.  Phys.  A  12,  619  (1979). 

[8]  Akira  Nakamura,  J.  Phys.  Soc.  Japan  47,  1335  (1979). 

[9]  T.  L.  Boch  and  M.  D.  Kruskal,  Phys.  Lett.  74A,  173  (1979). 

[10]  Y.  Matsuno,  J.  Phys.  A  13,  1519  (1980). 

[11]  H.  H.  Chen  and  D.  J.  Kaup,  Phys.  Fluids  23,  235  (1980). 

[12]  M.D.  Spector  and  T.  Miloh,  SIAM  J.  Appl.  Math  54,  688  (1994). 

[13]  Y.  Matsuno,  Phys.  Rev.  £?51,  1471  (1995). 

[14]  Y.  Matsuno,  Phys.  Rev.  E  52,  6333  (1995). 

[15]  M.J.  Ablowitz,  D.J.  Kaup,  A.C.  Newell  and  H.  Segur,  Stud.  Appl.  Math.  53,  249 
(1974). 

[16]  R.R.  Coiffman  and  M.V.  Wickerhauser,  Inverse  Problems  6,  825  (1990). 

[17]  R.L.  Anderson  and  E.  Tallin,  Lett.  Math.  Phys.  9,  299  (1985). 

[18]  N.  Levinson,  Danske  Vid.  Selsk.  Math.  Fys.  Medd.  25,  1  (1949).  See  also  G.L. 
Lamb,  Jr.,  ’’Elements  of  Soliton  Theory”,  p.  59  [John  Wiley,  New  York  (1980)]. 

[19]  P.  M.  Santini,  M.  J.  Ablowitz,  and  A.  S.  Fokas,  J.  Math.  Phys.  25,  892  (1984). 

[20]  See  Eq.(3.7)  in  Y.  Matsuno,  J.  Phys.  Soc.  Japan  51,  3734  (1982). 

[21]  A.S.  Fokas  and  B.  Fuchssteiner,  Phys.  Lett.  86 A,  341  (1981). 

[22]  J.D.  Meiss  and  N.R.  Perira,  Phys.  Fluids  20,  700  (1978). 

[23]  K.  M.  Case,  J.  Math.  Phys.  20,  972  (1979). 

[24]  H.  H.  Chen  and  D.  J.  Kaup,  J.  Math.  Phys.  21,  19  (1980). 

[25]  Y.  Matsuno,  Phys.  Lett.  87A,  15  (1981). 


22 


[26]  Y.  Matsuno,  J.  Phys.  Soc.  Japan  51,  667  (1982). 

[27]  Y.  Matsuno,  Intemat.  J.  Mod.  Phys.  B  9,  1985  (1995). 

[28]  Y.  Kodama,  M.  J.  Ablowitz,  and  J.  Satsuma,  J.  Math.  Phys.  23,  564  (1982). 


23 


